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Abstract 

We argue that the sharp-cutoff Wilson renormalization group provides a 
powerful tool for the analysis of second-order and weakly first-order phase 
transitions. In particular, in a computation no harder than the calculation of 
the 1-loop effective potential, we show that the Wilson RG yields the fixed point 
couplings and critical exponents of 3-dimensional 0(N) scalar field theory, with 
results close to those obtained in high-order e-expansion and large- N calculations. 
We discuss the prospects for an even more precise computation. 
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Weakly first-order and second-order phase transitions have been important for some 
time in cosmology. Particular attention has been focused on the electroweak phase 
transition (EWPT), in an effort to determine whether it is first-order and thereby capable 
of generating the observed asymmetry between matter and antimatter in the universe 
However, progress has been complicated by infrared problems that invalidate the usual 
loop expansion ||. In statistical mechanics, the usual remedy is to use the e-expansion 
[||, |], and this has been applied to the EWPT J5J. However, it has been argued that 
the e-expansion may be unreliable for the EWPT |J, since there may be a fixed point in 
three dimensions that is not visible in an expansion around four dimensions. In this letter 
we will show how a second-order phase transition can be analyzed by the 'Wilson" or 
"exact" renormalization group (RG) 0, || directly in three dimensions, without recourse 
to the e-expansion. It appears straightforward to extend this method to first-order phase 
transitions [I4[, which would enable us to determine the order and the dynamics of the 
EWPT and other cosmological phase transitions. 



There has recently been a resurgence of interest in the Wilson RG 
An important contribution was made by Tetradis and Wetterich [15|, who used the 
smooth-cutoff Wilson RG to calculate critical exponents of a second-order phase 
transition, and found impressive agreement with the traditional methods mentioned 
above. Unfortunately, the use of the smooth cutoff greatly complicates the computation, 
and numerical methods are required to determine the differential flow equations, which 
are then solved numerically to obtain the exponents. We will find, in contrast, that the 
sharp-cutoff Wilson RG yields the fixed point couplings and critical exponents in a simple 
and almost completely analytic form, requiring no more effort than the calculation of the 
one-loop effective potential. Even with a radically truncated action, the resultant values 
are within about 10% of those found by the traditional methods. In the appendix we 
discuss the problems that arise with both smooth and sharp cutoffs in a more precise 
calculation. 

The Wilson RG follows the flow of the effective action (free energy) as degrees of 
freedom are integrated out in successively shrinking momentum shells. This corresponds 
to calculating the effective action with an infra-red (IR) cutoff A for the loop corrections, 
and sending A — > 0. The effective action then flows with A, its couplings obeying 
differential equations that we will study below. For high temperature systems, we are 
interested in the flow for the 3 dimensional theory, and if there is a second-order phase 
transition we expect to find a fixed point that is IR-attractive in all coupling directions 
except one. This corresponds to the fact that one linear combination of the UV couplings 
(i.e. the temperature) must be fine-tuned to a critical value if the theory is to flow into 
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the fixed point in the IR. For temperatures close to the critical one, the effective action 
stays very close to the critical point for a large range of A, and this causes the singular 
dependence of many quantities on T — T c . 

We will consider the O(N) scalar field theory in 3 dimensions. It will turn out that 
on the critical surface the mass-squared flows to zero from below, so the true minimum is 
at a nonzero field. By the O(N) symmetry we can choose this to lie in the 4>i direction, 
so we shift the field (pi = <&5n + ipi, and the effective action for the shifted field ip is 

T = \zd^ t d^ t + v{ P ) 

p = = Po + $<Pl + \PMi 

This effective action is parameterized by the wavefunction renormalization Z, the 
position po = of the minimum of the potential, and the derivatives of V at p Q , 
V n = d n V/dp n (p ). By definition of p , Vl = 0. 

The procedure for obtaining the Wilson RG equations is simply to calculate the 
effective action in the usual way, but with an IR cutoff A in the propagators. We then 
differentiate with respect to A to obtain the flow of Z, po, and the V n with A [§, |11|, [12|, |I~3" 



The flow of Z is given by the anomalous dimension rj, which is discussed in the 
Appendix. For the sake of computational simplicity, we have assumed that Z is a 
constant, so 77 = 0. This is expected to be a good approximation, since e- and large-A?" 
expansions indicate rj ^ 0.04. The flow of the effective potential just corresponds to 
differentiating with respect to an IR cutoff A in the momentum integral of the usual 
1-loop formula flnj, so we easily obtain the flow equation, first written down by Nicoll, 
Chang, and Stanley (Ref. §, 1976): 



A dV(p) 



OA 



K 3 ^ 

-A 3 

2 



\n(ZA 2 + V\p) + 2pV"(p)) + (N — 1) ln(ZA 2 + V'(p) 



(2) 



where K 3 = l/(2n 2 ). We emphasize that although this looks like a 1-loop result there 
are no higher loop terms: the two-loop contribution from a given shell goes as (dA) 2 and 
does not affect the A-derivative. The approximation in Eq. (Q) is that it uses a truncated 
effective action, only taking account of the flow of the momentum-independent couplings 
(the effective potential). Including momentum-dependent vertices would clearly improve 
the accuracy of the results. 

Since we expect the dimensionless couplings to flow to a fixed point, we now scale 
out the dimensions and Z-dependence, defining a dimensionless field and couplings: 

p = p-, v(p) = A~ 3 V = v + Vl {p - po) + ±v 2 (p - p ) 2 + ■ ■ ■ (3) 
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Substituting these into Eq. @ we find the flow of the dimensionless couplings at constant 
p. Remembering that p is also flowing, our final result is 

A s = A 1 1 - *■+ % ' \ K H + + + A \ <«) 

We determine the flow of po by noting that by definition of po, v\ — for all A, so the 
flow equation for v\ becomes a flow equation for po ; 

a ^ = - 1a ^ - -» + - + 3 + (5) 

dA -u 2 dA 2 V 1 + 2p t> 2 / 

We now expand Eq. @ in powers of (p — p ) to find the flow of v 2 , v 3 , etc. 

A dv 2 . rfpo 

A— = v 3 A— -v 2 + pov 3 

l K, ((N - !)(«, - , 2 2 ) - + + ^ + 2P ° V4 



2 V /v /y (l + 2p ^) 2 l + 2p v 2 

dv 3 dp _ 

A— = , 4 A— + p , 4 (6) 



2 °V Jy z *> (l + 2p v 2 f 

(2>v 2 + 2p Q v 3 )(5v 3 + 2p v 4 ) 7v 4 + 2p v 5 
(1 + 2p ^) 2 (1 + 2p v 2 ) 

Expanding to arbitrary order to find the flow of v n is a matter of tedious but 
straightforward algebra. 

To find the fixed point we set all flows to zero and solve Eqs. fl5|) and for the fixed 
point couplings Pq,v 2 ,v 3 . . .v*. If our results are to be trustworthy then the fixed point 
couplings must settle down at constant values as we relax the truncation, i.e. as n — > oo. 
Fortunately, finding the fixed point always involves solving two simultaneous equations, 
no matter what the value of n, since the po equation eliminates v 3 in terms of pQ and v 2 , 
the v 2 equation eliminates v\, and so on. We finally set v * +1 and v * +2 to zero, obtaining 
two simultaneous (polynomial) equations for pQ and v 2 . Solving these is the only part of 
the calculation which we have to perform numerically. The results for the two scalar field 
case (N = 2) are shown in Table 1 as a function of n. (They were obtained using the 
symbolic manipulation program Mathematica to analytically expand the flow equations 
and numerically find pQ and v 2 .) We see impressive convergence of the couplings and the 
critical exponent (3 (calculated below) as the truncation of the potential is relaxed. 

To calculate critical exponents, consider what happens when T is not quite equal to 
T c , so that the dimensionless couplings almost flow into the fixed point, but turn away at 
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Table 1: The N = 2 fixed point couplings and critical exponent (5, calculated using n 
evolution equations, i.e. including vertices up to <p 2n . Note the striking convergence with 
increasing n. (M = 10 6 ). 



the last moment and explode exponentially away from it. From Eq. (HD we see that when 
Po becomes of order 1, at some A = A/, the noncanonical scaling term (multiplied by Ks) 
becomes negligible, and from then on po grows canonically, so po tends to a constant: 



Po 



A, 



(7) 



The critical exponent (3, defined by p ~ (T — T c ) 2/3 , can then be found by examining 
the flows in the neighborhood of the fixed point. Write the couplings as a vector v, with 
components Vi — v * (let V\ be po for this purpose); then the flows can be linearized: 



A- 



dv 



W - v 



so that 



W, 



_d_ A d(vi-v*) 



dA 



(8) 



(9) 



We will find that W has one positive eigenvalue wi, corresponding to the repulsive 
direction, with eigenvector wi. It follows that 



v(A) = A(T - T C )A- Wl w 1 + ■ ■ 



(10) 



The couplings leave the vicinity of the fixed point when the coefficient of Wi becomes of 



order 1, so A(T — T C )A 



— Wl 

f 



' 1. By Eq. (0), we otain the critical exponent (3: 
(T_ Tc )iM so p = l/(2wi), 



(11) 
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Table 2: The O(N) fixed point couplings and critical exponents for various N. The 
values in square brackets are the current best estimates, based on e-expansion to order 
e 5 (N < 3) ||, or large N to order 1/iV 2 (N > 10) [0. 



where W\ is the largest eigenvalue of the linearized flow equations ([5]) at the fixed point. 
We can now calculate the fixed point couplings, and the critical exponent f3 for any N . 
(Except for the anomalous dimension, the other exponents follow by scaling relations.) 
The results are given in Table 2. 

We see that there is striking agreement between the sharp-cutoff Wilson RG, and 
conventional e-expansion and large- iV results, particularly so considering the crudity of 
our approximation, in which all momentum dependence in the effective action has been 
ignored. These results also agree closely with those obtained numerically by Tetradis and 
Wetterich |15] via the smooth-cutoff Wilson RG. It seems reasonable to expect that the 
inclusion of momentum- dependent vertices would enhance the accuracy of the results still 
further. However, as we show in the appendix, this is not a completely straightforward 
procedure. In particular, the smooth-cutoff calculation |T5|] of the anomalous dimension 
appears to depend on the profile of the cutoff, rendering its physical significance unclear. 

For cosmological applications, it is necessary to formulate a gauge-invariant version 
of the momentum cutoff, in order to study gauge theories with weakly first-order phase 
transitions such as the electroweak phase transition. Progress in this direction has already 
been reported |]18|| . It would also be interesting to see if the sharp-cutoff Wilson RG 
could be applied with equal success to second order phase transitions in other statistical 
mechanical models, or to calculation of the effective potential in other field theories where 
perturbation theory breaks down in the IR, such as multiple scalar field theories with 
radiatively induced first-order phase transitions of the Coleman- Weinberg type [14]. 
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Figure 1: The diagram that contributes to the coefficient of k 2 in the 2-point 1PI Green's 
function, i.e. to the running of Z. 
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Appendix A The anomalous dimension 

The anomalous dimension is given by the flow of the wavefunction renormalization Z, so 
it is convenient to express the effective action in terms of ipi rather than p: 

v(p) = ($V? + &<pi<pm + l<Pi<pm<Pj) + ■■■ ( A - x ) 

The first term is the mass of the radial mode, which comes into our flow formulae in 
its dimensionless form, as 2p t>2- The next term provides a cubic interaction between 
the Goldstone and radial modes. This is the term that, via the graph shown in Fig. 1, 
will drive the cutoff-dependence of the wavefunction renormalization Z of the Goldstone 
modes. 

We calculate Z = ^d^^/dk 2 by evaluating this diagram in the usual way, but with 
an IR cutoff 6 e (p 2 — A 2 ) (which becomes a step function as e — > 0) with every propagator. 
Differentiating with respect to A, 



= ^ = ( ^v) 2 --— f d3q I ^(g~ A ) e £ {{ q + k) 2 -K 2 ) 

% ZdA 1 2) Z2dk 2 J (2tt) 3 \Zq 2 + 2p V 2 Z(q + k) 2 

A) 6 £ {{q + k) 2 -\ 2 ) \ 
Z(q + k) 2 + 2p V 2 j 




(A.2) 



This formula is to be evaluated in the limit k ^ 0, and the sharp cutoff limit e — ► 0; 
but in which order? Clearly if we send k — > first then the fc-derivatives convert the 
^-functions into <5 e -functions (and their derivatives) centered at q — A, and we obtain a 
result that is sensitively dependent on e, the profile of the cutoff, and diverges as e — > 0. 
Moreover, the expansion in k is only valid for k ~ e |R|. This is not a problem with the 
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sharp cutoff limit, as has been suggested [|i~0|, 11|, since it is just as bad for the smooth 



cutoff to give profile- dependent results as it is for the sharp cutoff to give infinite ones. 
Rather, it implies that the order of limits is wrong. Physically, we want our effective 
action to be valid for momenta spanning as much of the range to A as possible. This is 
achieved by sending e — > first, and we then obtain e-independent results for any k 3> e, 
just as one would expect. 

A more serious problem is that the resultant 2-point function contains a non-local 
\k\ term, as well as the expected k 2 term (see also Ref. [H|). This occurs because we 



have the same infrared cutoff on two different loop momenta, and only part of the loop 
is integrated out. This does not mean that the Greens functions, which are physically 
measurable, are non-local. The effective action with an IR cutoff is not a physically 
measurable quantity, and so it may indeed contain non-local terms It is not clear 
whether we should include a \k\ term in our effective action, with a separate flow equation 
for its coefficient, or whether we can amend the definition of the coarse-grained effective 
action in such a way as to eliminate non-local terms. This issue clearly deserves further 
study. 
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